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2 套 气象 数据 在 资料 缺乏 地 区 的 适用 性 评估 
一 一 以 呼 图 壁 河流 域 为 例 
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摘 要 : 在 无 资料 地 区 ,降水 .气温 等 观测 资料 缺乏 ,影响 洪水 预报 的 准确 性 ,进而 影响 防汛 抗旱 ,水 资源 规划 等 工 
作 的 开展 。 因 此 ,研究 无 资料 地 区 已 有 降水 气温 数据 的 适用 性 分 析 十 分 必要 。 应 用 包含 融雪 径流 模拟 的 HBV 模 
型 ,在 呼 图 辟 河 流域 石门 水 文 站 以 上 区 域 ,基于 模拟 出 的 融雪 洪水 径流 与 实测 径流 对 比 ,分 别 对 中 国 地 面 降 水 与 气 
温 日 值 0.5°x0.5° 格 点 数据 集 和 气象 站 观测 数据 集 进行 了 对 比分 析 及 其 融雪 型 径流 模拟 的 适用 性 分 析 。 结 果 表 明 : 
格 点 数据 驱动 下 水 文 模型 的 融雪 径流 模拟 效果 总 体 优 于 利用 站 点 观测 数据 驱动 的 水 文 模型 。 格 点 数据 驱动 下 水 
文 模型 模拟 融雪 径流 的 Nash 系数 在 验证 期 为 0.792, 站 点 数据 驱动 下 水 文 模型 模拟 融雪 径流 的 Nash 系数 在 验证 期 
为 0.433。 同 时 分 析 了 呼 图 壁 河 流域 融雪 洪水 的 一 些 特征 ,并 结合 驱动 数据 的 不 同 ,分 析 了 模拟 融雪 径流 出 现 误 差 
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的 可 能 原因 ,为 如 何 提高 具有 融雪 特征 的 无 资料 地 区 融雪 型 洪水 预报 准确 性 提供 支撑 。 


关键 词 : 降水 数据 ; 无 资料 地 区 ; 融雪 洪水 ; 径流 模拟 ; 适用 性 评估 ; 呼 图 壁 河流 域 


在 融雪 融 冰 地 区 ,进行 准确 的 径流 模拟 与 预报 
对 水 资源 规划 与 管理 ,减少 融雪 洪水 带 来 的 损失 十 
分 重要 。 大 多 数 融 雪 过 程 发 生 在 高 寒 地 区 或 河 
流 源头 ,这 些 地 方 往往 缺乏 实测 资料 ”。 因 此 ,对 于 


冰川 与 积 雪 ,融化 后 会 产生 融雪 洪水 ""。 近 20 a 来 ， 
新 疆 融 雪 ( 冰 ) 型 洪水 占 新 疆 洪水 的 39%"”。 现 阶 
段 ,我 国 的 国家 气象 站 基本 保持 在 2200~2400 个 左 
^h ,东南 地 区 分 布 较 多 ,西北 地 区 分 布 较 少 呈 。 新 


这 些 具 有 融雪 特征 的 资料 缺乏 地 区 ,研究 如 何 提高 
水 文 预报 准确 性 十 分 必要 。 在 高 寒山 区 ,降水 特征 
受 地 形 影 响 较 大 ,复杂 的 山区 地 形 和 局 地 小 气候 形 
成 了 这 异 多 变 的 降水 时 空 分 布 格局 “。 气 温 的 分 布 
相当 程度 上 也 受到 地 形 因 素 的 影响 ,具有 海拔 引起 
的 垂直 差异 ”。 在 寒 区 水 文 模拟 过 程 中 ,降水 和 和 气 
温 是 关键 的 变量 数据 "。 在 融雪 洪水 模拟 中 ,降水 


疆 地 区 的 气象 站 总 体 较 为 稀 梳 。 在 呼 图 壁 河流 域 
石门 水 文 站 以 上 的 区 域 ,全 国 2400 个 国家 级 地 面 气 
象 观 测 站 中 没有 站 点 在 该 区 域 中 。 除 气象 站 点 数 
据 外 ,多 机 构 发 布 了 多 种 不 同 的 降水 、 气 温 数 据 产 
品 ,不 同 数据 产品 在 数据 来 源 、. 时 间 序 列 长 度 .数据 
处 理 方式 与 数据 用 途上 各 有 不 同 “。 对 资料 缺乏 
地 区 ,研究 已 有 降水 ` 气温 数据 的 适用 性 十 分 必 


和 气温 数据 是 影响 水 文 模型 输入 不 确定 性 的 重要 
ust. Bee 7K BE ae SERI CARCER ,往往 是 水 文 模 
型 模拟 径流 时 最 重要 的 输入 数据 。 因 此 ,对 不 
同 来 源 的 驱动 数据 ,如 降水 、 气 温 数 据 进 行 适 用 性 
评估 ,对 提高 预报 的 准确 性 有 重要 意义 。 

新 疆 位 于 中 国 西 北 ,是 典型 的 内 陆 干 旱 区 ,内 
有 许多 高 大 山体 ,在 低温 条 件 下 ,高 山区 降水 形成 
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要 。 国 家 气象 中 心 提 供 的 中 国 地 面 气温 与 降水 日 
值 0.5°x0.5° 格 点 数据 集 包 括 1961 年 至 今 的 逐日 气 
温 与 降水 数据 ,时 间 序 列 长 度 足 够 , 且 该 数据 是 基 
于 国家 级 高 密度 台 站 插值 而 来 ,是 中 国 经 纬度 边界 
范围 内 的 气温 、 降 水 数据 , 比 诸多 全 球 范 围 的 气温 、 
降水 产品 更 有 针对 性 的 。 网 格 降水 数据 产品 可 有 
效 地 替代 地 面 实测 数据 ,但 需要 经 过 评估 验证 5 。 
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Gao 等 “以 基于 站 点 的 中 国 格 网 日 降水 分 析 产 品 
CGDPA 为 基准 ,评估 了 2 套 降 水 数据 产品 在 2 处 无 
资料 地 区 的 精度 ,并 驱动 水 文 模型 分 析 了 其 对 洪水 
频率 分 析 的 适用 性 。Dembélé 等 "在 评估 总 结 多 套 
不 同 降水 和 气温 数据 产品 的 组 合 驱 动 mHM 水 文 模 
型 的 效果 时 ,发 现 不 存在 始终 保持 最 优 的 气温 或 降 
水 数据 集 ,一 些 区 域 性 数据 产品 优 于 全 球 数 据 产 
Mo Mazzoleni 等 在 研究 不 同 降水 数据 集 驱 动 分 
布 式 水 文 模型 模拟 径流 时 发 现 , 不 同 的 流域 有 不 同 
的 最 优 降 水 数据 ,模拟 结果 对 流域 特征 十 分 敏感 。 
在 以 往 人 研究 中 ,研究 区 域 的 空间 范围 也 往往 较 大 。 
因此 ,对 于 具体 的 中 小 无 资料 流域 ,针对 各 自 不 同 
的 情况 , 仍 需要 因地制宜 ,做 进一步 的 分 析 与 评 佑 。 

HBV (Hydrologiska Byráns Vattenbalansavdeln- 
ing) 模 型 是 瑞典 气象 和 水 文 研究 所 (SMHIT) 于 20 世 
纪 70 年 代 开 发 的 概念 性 水 文 模型 ”。 对 融雪 径流 
作用 机 制 明显 的 流域 ,HBV 模 型 的 适应 性 较 强 。 
目前 ,HBV 模 型 已 被 应 用 于 全 球 50 多 个 国家 ,被 广 
泛 应 用 于 各 种 尺度 ,成 为 许多 应 用 类 型 中 更 普遍 
的 工具 ,在 许多 无 资料 或 资料 缺乏 的 地 区 也 得 到 
nga. 

本 文 以 气象 站 点 稀疏 的 新 疆 呼 图 壁 河流 域 厂 
门 站 以 上 的 区 域 为 研究 区 域 ,选取 中 国 气象 局 的 气 


象 站 观测 数据 集 和 中 国 地 面 降水 与 气温 日 值 0.5°x 
0.5" 格 点 数据 集 ,应 用 包含 融雪 径流 模拟 的 HBV FR 
型 ,在 呼 图 壁 河流 域 石门 水 文 站 以 上 区 域 , 基 于 模 
拟 出 的 融雪 洪水 径流 与 实测 径流 对 比 ,分 别 对 选取 
的 2 套 气 温 、 降 水 数据 集 进 行 对 比分 析 及 其 融雪 型 
径流 模拟 的 适用 性 分 析 , 分 析 了 呼 图 壁 河流 域 模拟 
与 实测 融雪 洪水 的 一 些 特征 ,也 结合 驱动 数据 的 不 
同 分 析 了 模拟 融雪 径流 出 现 误 差 的 可 能 原因 ,为 如 
何 提高 无 资料 地 区 融雪 型 洪水 预报 的 准确 性 提供 
文 撑 。 


1 数据 和 方法 


1.1 研究 区 概况 

本 文人 研究 区 域 为 新 疆 呼 图 壁 河流 域 石门 站 以 
上 的 区 域 (图 1), 呼 图 壁 河流 域 位 于 天 山北 坡 中 段 、 
准噶尔 盆地 南 缘 ,地理 位 置 介 于 86*05'~87"08' 下 ， 
43°07'~45°20'N 之 间 。 流 域 自 上 而 下 有 石门 水 文 
站 、 青 年 渠 首 及 茂 茂 坝 分 水 曾 站 。 石 门 站 以 上 流域 
集 水 面积 为 1840 km’”, 河 道 坡 降 达 23.3%o, 年 径流 量 
占 呼 图 壁 河 全 流域 年 径流 量 的 93.3%"。 呼 图 壁 河 
径流 年 内 分 配 极 不 均匀 ,在 夏季 相对 集中 , 占 全 年 
径流 量 的 80%。 流 域 地 表 径 流 补 给 主要 来 源 于 高 
山 雪 冰 融 水 、 季 节 性 积 雪 融 水 降水 和 地 下 水 *。 
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注 :图 中 经 纬度 网 格 线 亦 为 格 点 单元 分 割 线 。 


图 1 研究 区 示意 图 


Fig.1 Schematic diagram of the study area 
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1.2 HBV 模 型 

HBV 模型 结构 简单 ,易于 理解 。 本 文 使 用 的 
HBV 模型 版 本 为 AghaKouchak 等 "使 用 的 HBV 集 
总 式 日 模型 版 本 。 值 得 注意 的 是 ,因为 研究 区 域 气 
象 站 点 稀 玖 ,没有 足够 的 蒸 散 发 数据 作为 驱动 数 
据 ,因此 ,在 求 潜在 蒸 散发 时 ,将 HBV 模 型 中 原本 用 
月 平均 气温 、 月 平均 潜 蒸 散发 .日 平均 气温 求 修正 
的 潜 蒸 散发 的 方法 改 为 使 用 修正 的 Blaney-Criddle 
方程 ,根据 风速 .湿度 .太阳 辐射 等 因素 ,以 日 平均 
气温 为 主要 驱动 数据 计算 修正 的 潜在 蒸 散发 量 。 
修改 后 ,模型 输入 的 驱动 数据 为 日 降水 量 和 日 平均 
气温 数据 ,模型 的 输出 为 日 平均 流量 。 

模型 总 体 水 平衡 可 描述 为 : 


P-E-R=§(SP+SM+ UZ-LZ4L,.) (1) 


式 中 :P 是 降水 量 ; 到 是 蒸 散 发 量 ;及 是 径流 量 ;SP 为 

REA mE; SM 为 土壤 含水 量 ; UZ 为 表层 地 下 含水 
Ht LZ 为 深层 地 下 含水 量 ; Li 为 流域 中 水 体 的 
水 量 。 

HBV 模型 包含 4 个 主要 模块 :融雪 与 积 雪 模块 ; 
土壤 水 分 和 有 效 降 水 模块 ; 蒸 散发 模块 ;径流 响应 
模块 。 

(1) 融雪 与 积 雪 模 块 

HBV 模型 假定 融雪 和 积 雪 由 观测 的 气温 (7) 和 
闷 值 温度 (77) 之 间 的 相对 关系 控制 。 如 果 7T<77,， 
则 假定 所 有 进入 的 降水 都 以 降雪 的 形式 积累 ,不 产 
生 径 流 。 如 果 了 >77, 则 积 雪 会 融化 ,融雪 用 度 日 法 
计算 。 液 态 水 量 为 降水 量 与 融雪 所 得 水 量 之 和 ,并 
随后 划分 为 入 渗 和 径流 成 分 。 度 日 法 计算 公式 为 : 

S, - DD(T- TT) (2) 
CP S, NME Emm: d) ; DD 2 BE H PSI f (mm Ct 
d;7 为 日 平均 气温 (%C);77 为 阔 值 温度 (%D)。 

(2) 土壤 水 分 和 有 效 降 水 模块 

根据 降水 时 土壤 含水 量 估算 ,模型 将 降水 分 为 
2 个 部 分 ,一 部 分 有 助 于 渗透 到 土壤 中 , 男 一 部 分 有 
助 于 产生 地 表 径 流 , 即 有 效 降 水 Pu。 


Bs (Mers (3) 


式 中 :Pr 表 示 有 效 降水 (mm);SM 是 实际 土壤 合 水 量 
(mm); FC 是 最 大 土壤 持 水 量 (mm);P 是 日 降水 量 
(mm);B 是 形状 因子 。 

(3) 蒸 散发 模块 


实际 蒸 散 发 量 由 实际 土壤 含水 量 与 PWP 的 关 
系 经 以 下 的 计算 得 到 : 


T Ip) SM « PWP 

Hs , SM» PWP 

式 中 :zw 表示 潜在 蒸 散发 nm); 及 表示 实际 的 蒸 散 

发 量 (mm); PP 表示 实际 蒸 散 发 达到 潜在 蒸 散 发 时 
的 土壤 含水 量 (mm)。 

(4) 径流 响应 模块 

HBV 模型 利用 线性 水 库 的 概念 来 模拟 估算 流 
域 出 口 处 的 径流 。 系 统 由 2 个 虚拟 的 水 库 组 成 ,第 1 
个 水 库 用 来 描述 近 地 表 径流 过 程 ,第 2 个 水 库 用 来 
模拟 基 流 。 这 2 个 水 库 通过 恒定 的 渗透 速率 相互 连 
接 。 总 的 模拟 径流 量 通过 将 第 1 个 和 第 2 个 水 库 流 
出 量 相 加 得 出 。 

1.3 数据 来 源 

采用 的 数据 集 为 中 国 气象 局 的 气象 站 公布 的 2 
套数 据 集 。 第 1 套数 据 为 站 点 观测 数据 ,选取 全 国 
2400 国家 级 地 面 气象 观测 站 中 的 3 个 站 点 和 石门 水 
文 站 的 观测 数据 。 在 全 国 2400 国家 级 地 面 气象 观 
测 站 中 没有 气象 站 在 研究 区 域内 ,选取 距离 研究 区 
域 最 近 的 天 山大 西 沟 站 与 较 近 的 呼 图 壁 站 、 小 汇 子 
站 ,分 别 距 研究 区 边界 17 km、39 km、45 km 左右 (图 
1)。 对 选取 的 4 个 站 点 求 其 气温 、 降 水 数据 的 算术 
平均 值 ,作为 第 1 套 驱 动 数 据 。 

第 2 套数 据 为 气温 与 降水 格 点 数据 ,降水 格 点 
数据 来 自 中 国 地 面 降水 日 值 0.5°x0.5° 格 点 数据 集 
(V2.0) ,数据 集 代码 :SURF_CLI_CHN_PRE_DAY_ 
GRID_0.5。 气 温 格 点 数据 来 自 中 国 地 面 气 温 日 值 
0.5? x0.5° 格 点 数据 集 (V2.0) ,数据 集 代 码 : 
SURF_CLI_CHN_TEM_DAY_CRID_0.5。 选 取 本 文 
人 研究 区 所 在 的 3 个 格 点 单元 (图 1) ,并 按 流域 在 格 点 
单元 中 所 占 的 面积 比例 大 小 对 所 选 格 点 单元 的 降 
水 .气温 数据 加 权 平 均 ,将 得 到 的 降水 .气温 日 数据 
作为 男 1 套 驱动 数据 。 

本 文 的 研究 区 内 有 石门 水 文 站 ,位 于 人 研究 区 出 
口 。 石 门 水 文 站 的 气温 、 降 水 与 径流 数据 均 来 自 水 
文 年 鉴 。 

1.4 统计 指标 

DLH HBV 模型 对 研究 区 进行 融雪 洪水 模拟 , 率 
定期 为 2000 年 1 月 1 日 至 2007 年 12 月 31 日 ,验证 
期 为 2008 年 1 月 1 日 至 2011 年 12 月 31 日 。 


(4) 
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径流 模拟 的 结果 用 Nash 效率 系数 与 流量 平均 
相对 误差 判定 : 
240470.) 
28470 
m 


obs 
式 中 :0 为 实测 流量 s) ; Qus Ar RF Tit et [EL (m + 
ss Qu, ADULT Tie EF [B (m? +s") 5 0O 为 模拟 流 
EF X fi (m^ - s), NSE 29 Nash 效率 系数 (Nash- 
Suttcliffe Efficiency) , MRE 为 平均 相对 误差 (Mean 
Relative Error), Nash 系数 NSE 介 于 -co~1 之 间 , 其 
值 越 接近 1, 则 模拟 结果 越 精 确 。 流 量 平均 相对 误 
差 MRE 用 来 评估 总 体 流量 是 否 平衡 。 


2 结果 与 分 析 


2.1 站 点 数据 与 格 点 数据 的 对 比分 析 

首先 ,将 2 套 驱 动 数据 一 一 站 点 气温 .降水 数据 
与 格 点 气温 、 降 水 数据 做 对 比 (图 2)。 

从 图 2a 可 以 看 出 ,站 点 日 平均 气温 整体 高 于 格 


NSE=1- (5) 


MRE = x 100% (6) 


| ^ | 
Mi A i i i Il 


气温 /°C 
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点 日 平均 气温 。 从 2000 年 1 月 1 日 至 2011 年 12 月 
31 日 的 日 平均 气温 特征 值 来 看 ,2000 一 2011 年 ,站 
点 日 平均 气温 的 最 大 值 为 24.3 % ,在 2004 年 7 月 14 
H , 格 点 日 平均 气温 的 最 大 值 为 17.5% ,在 2006 年 7 
月 31 日 。 站 点 日 平均 气温 的 最 小 值 为 2011 年 1 月 7 
日 的 -24.8 C , 格 点 日 平均 气温 的 最 小 值 为 2006 年 1 
月 4 日 的 -26.6 %C。 

从 季节 尺度 看 ,在 2000 年 1 月 1 日 至 2011 年 12 
月 31 日 ,夏季 站 点 日 平均 气温 数据 的 多 年 平均 值 为 
16.1% , 比 格 点 日 平均 气温 数据 的 夏季 多 年 平均 值 
9.7 避 高 6.4% ,在 4 个 季节 中 相差 最 大 。2 套 日 平均 
气温 数据 的 春季 多 年 平均 值 之 差 次 之 ,春季 站 点 日 
平均 气温 数据 的 多 年 平均 值 为 44 %C, 比 格 点 数据 
的 -0.8 CC 高 5.2 %。 秋 季 的 情况 与 春季 较为 接近 , 秋 
季 站 点 日 平均 气温 的 多 年 平均 值 为 3.8 %C, 比 格 点 
数据 的 -0.6 % 高 44%。 在 4 个 季节 中 ,冬季 2 套 日 
平均 气温 数据 的 多 年 平均 值 相差 最 小 ,站 点 日 平均 
气温 数据 的 冬季 多 年 平均 值 为 -11.7 %C, 比 格 点 日 平 
均 气 温 数据 的 -13.4% 高 1.7%C。 

从 多 年 尺度 看 ,在 2000 年 1 月 1 日 至 2011 年 12 


一 一 气温 ( 格 点 数据 ) 一 一 气温 (站 点 数据 ) 
I! Ki I "i 
| | ' lj i) j i " | 


1 | i m | | | 
a NA. | | 
4 Į) 


! 
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降水 量 /mm 
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一 一 降水 量 ( 格 点 数据 ) 
—— 降水 量 (站 点 数据 ) 
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图 2 站 点 数据 与 格 点 数据 对 比 


Fig.2 Comparison of station data and grid data 
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月 31 日 ,站 点 日 平均 气温 数据 的 多 年 平均 值 为 
3.2 % , 格 点 日 平均 气温 数据 的 多 年 平均 值 为 
-1.2 % ,站 点 日 平均 气温 数据 的 多 年 平均 水 平 比 格 
点 数据 高 4.4 C。 

分 析 图 2b 可 知 ,在 每 年 日 降水 量 最 大 的 几 天 ,2 
种 降水 数据 产品 有 较 大 偏差 ,所 用 站 点 的 日 降水 数 
据 往往 小 于 格 点 日 降水 数据 。 以 2000 年 为 例 , 在 
2000 年 ,站 点 数据 和 格 点 数据 的 日 降水 量 最 大 值 均 
在 6 月 11 日 ,这 一 天 站 点 数据 日 降水 量 为 18.7 mm, 
格 点 数据 日 降水 量 为 39.9 mm。 在 2000 年 ,站 点 数 
据 最 大 3 日 降水 量 为 36.5 mm, 格 点 数据 最 大 3 日 降 
水 量 为 65.0 mm。 站 点 数据 中 ,根据 24h 降 水 量 , 属 
于 中 雨 (10~25 mm) 的 有 9 d, 属 于 小 雨 (0.1~10 mm) 
的 有 171 d, 其余 (<0.1 mm) 的 有 186 d。 格 点 数据 
中 ,根据 24b 降 水 量 , 属 于 大 十 (25~50 mm) 的 有 
2 d, J&& T "P RR (10-25 mm) 的 有 13 d, 属 于 小 雨 (0.1~ 
10 mm) 的 有 234 d, 其 余 (<0.1 mm) 的 有 117 do 

从 季节 尺度 来 看 ,在 春季 ,站 点 数据 日 降水 量 
平均 值 为 1.2 mm, 格 点 数据 日 降水 量 平 均值 为 1.0 
mm。 在 夏季 ,站 点 数据 日 降水 量 平均 值 为 2.4 mm, 
格 点 数据 日 降水 量 平均 值 为 3.5 mm。 在 秋季 ,站 点 
数据 与 格 点 数据 日 降水 量 平均 值 均 为 0.7 mm. 在 
冬季 ,站 点 数据 日 降水 量 平均 值 为 0.3 mm, 格 点 数 
据 日 降水 量 平均 值 为 0.2 mm。 

从 多 年 尺度 来 看 ,在 2000 年 1 月 1 日 至 2011 年 
12 月 31 日 ,站 点 数据 的 多 年 平均 降水 量 为 426.3 
mm, 格 点 数据 的 多 年 平均 降水 量 为 494.2 mm。 格 
点 降水 数据 总 体 大 于 站 点 降水 数据 。 

从 2000 年 1 月 1 日 至 2011 年 12 月 31 日 各 种 降 
水 类 型 的 天 数 看 , 格 点 数据 日 降水 量 最 大 值 为 39.9 
mm, TR Hi 24 了 降 水量, 属于 大 雨 (25~50 mm) 的 有 
14 d, 属 于 中 雨 (10~25 mm) 的 有 95 d, 属 于 小 雨 
(0.1~10 mm) 的 有 2771 d, 其余 (<0.1 mm) £f 1503 
d。 站 点 数据 日 降水 量 最 大 值 为 33.0 mm ,根据 24 h 
降水 量 ,属于 大 两 (25~50 mm) 的 有 1d, 属 于 中 雨 
(10-25 mm) 的 有 81 d, 属 于 小 雨 (0.1~10 mm) WA 
2067 d ,其余 (<0.1 mm) 有 2234 d, 

从 降水 的 空间 分 布 来 看 ,对 于 站 点 数据 ,2000 
年 1 月 1 日 至 2011 年 12 月 31 日 , 呼 图 壁 站 多 年 平均 
降水 量 为 210.0 mm ,小 渠 子 站 多 年 平均 降水 量 为 
586.5 mm, 天 山大 西 沟 站 多 年 平均 日 降水 量 为 
502.9 mm ,石门 站 多 年 平均 降水 量 为 405.8 mm。 如 


图 1, 结合 海拔 高 度 可 以 看 出 ,在 本 文 研 究 区 内 , 海 
拔 越 高 ,站 点 数据 的 多 年 平均 降水 量 越 大 。 对 于 格 
点 数据 ,2000 年 1 月 1 日 至 2011 年 12 月 31 日 , 按 流 
域 在 格 点 单元 中 所 占 的 面积 比例 大 小 加 权 之 前 , 左 
上 格 点 单元 (43.5"~44*N , 86? -86.5* E) I] EY 
降水 量 为 542.2 mm, 左 下 格 点 单元 (43°~43.5°N， 
86°~86.5°E) 的 多 年 平均 降水 量 为 469.0 mm ,右上 格 
点 单元 (43.5°~44°N,86.5°~87E) 的 多 年 平均 降水 
量 为 526.2 mm。 在 本 文 研 究 区 内 , 格 点 数据 的 多 年 
平均 降水 量 随 格 点 单元 的 平均 海拔 高 度 的 变化 规 
律 不 明显 。 

2.2 融雪 洪水 径流 模拟 

首先 ,采用 算术 平均 法 将 2000 年 1 月 1 日 至 
2007 ^F. 12 H 31 日 的 站 点 气温 、 降 水 数据 和 格 点 气 
温 、 降 水 数据 进行 融合 ,得 到 一 套 融 合 的 气温 .降水 
数据 ,用 以 率 定 模型 参数 。 

本 文 所 用 的 HBV 模型 共有 参数 10 个 ,其 中 度 
AAI DD 田间 持 水 量 FC HEMP YE 7k RB ks FEAL 
IRB, AR A URE BU, H FETS EAT A 
工 率 定 , 通 过 人 工 率 定 的 过 程 理解 每 个 参数 对 模拟 
结果 的 影响 ,再 采用 遗传 算法 对 敏感 参数 进行 率 定 
( 表 1)。 

使 用 率 定 好 的 参数 ,分别 用 站 点 数据 和 格 点 数 
据 驱 动 HBV 模型 模拟 融雪 洪水 的 结果 如 表 2 所 示 。 

分 别 绘制 融合 数据 与 站 点 数据 、 格 点 数据 驱动 
HBV 模型 模拟 的 融雪 洪水 流量 与 实测 流量 过 程 线 ， 
如 图 3 所 示 。 

使 用 融合 的 气温 .降水 数据 驱动 HBV 模型 率 定 
参数 , 率 定 期 融雪 洪水 模拟 的 Nash 系数 为 0.680, 平 
均 流 量 相 对 误差 为 1.63% , 率 定期 模拟 的 洪水 过 程 
如 图 3a 所 示 。 可 以 看 出 ,融合 数据 模拟 出 的 夏季 洪 
水 的 洪峰 部 分 普遍 偏 小 ,2001 年 .2003 年 .2004 年 、 
2005 年 .2006 年 .2007 年 的 春季 模拟 出 的 融雪 径流 
偏 大 ,会 模拟 出 实测 不 存在 的 洪水 过 程 。 在 部 分 年 
份 最 后 1 个 洪水 之 后 的 部 分 时 段 ,如 2000 年 9 月 至 
10 月 上 旬 、2001 年 10 月 中 旬 、2002 年 10 月 、2003 年 
10 H ,2004 4F 10 H ,2005 4F 10 H ,2007 4F 10 H , fat 
合 数据 驱动 HBV 模型 模拟 出 的 流量 减少 的 过 程 与 
实测 过 程 拟 合 较 好 。 在 流量 较 低 的 冬季 , 2004 年 
1 一 2 月 与 2005 年 2 月 融合 数据 驱动 HBV 模型 模拟 
出 的 流量 与 实测 流量 拟 合 较 好 ,其 余年 份 冬季 融合 
数据 驱动 HBV 模型 模拟 出 的 流量 普遍 低 于 实测 
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表 1 HBYV 模型 的 参数 范围 与 率 定 结果 


Tab.1 Parameters range and calibration results of HBV model 


模块 参数 参数 含义 参数 范围 率 定 结果 参考 文献 
融雪 与 积 雪 模块 TT [2d (Fil BE/°C -1.0-2.7 2.40 25] 
DD BE H ES (mmc. d^) 0.5-4.0 0.5 22] 
有 效 降 水 与 土壤 含 FC 最 大 土壤 持 水 量 /mm 100~450 132.88 25] 
水 量 模块 B 形状 系数 1-5 2.3 25] 
蒸 散发 模块 PP 蒸 散发 达到 潜在 蒸 散发 时 的 土壤 水 90-180 106 21] 
量 /mm 
径流 响应 模块 ko 地 表 径 流 退 水 系数 /d" 0.05~0.2 0.0500 21] 
lo 阔 值 水 位 /mm 2-5 2 21] 
kı 壤 中 流 退 水 系数 /d 0.01~0.2 0.0143 25] 
k 基 流 退 水 系数 /d" 0.001~0.1 0.0050 25] 
k, 渗透 系数 /d" 0.01~0.05 0.0200 21] 


表 2 使 用 不 同 驱动 数据 的 流量 模拟 结果 
Tab.2 Flow simulation results using 


different driving data 


时 期 驱动 数据 年 份 NSE MRE/% 
率 定期 融合 数据 ”2000 一 2007 0.680 1.63 
验证 期 站 点 数据 ”2008 一 2011 0.433 14.64 
格 点 数据 ”2008 一 2011 0.792 2.44 
y LE o 


用 站 点 数据 驱动 HBV 模型 模拟 融雪 洪水 过 程 
(图 3p) ,验证 期 融雪 洪水 模拟 的 Nash 系数 为 0.433， 
平均 流量 相对 误差 为 14.64%。 用 格 点 数据 驱动 
HBV 模型 模拟 融雪 洪水 过 程 (图 3c) ,验证 期 融雪 洪 
水 模拟 的 Nash 系数 为 0.792 ,平均 流量 相对 误差 为 
2.44%。 从 Nash 系数 与 平均 相对 误差 来 看 , 格 点 数 
据 驱 动 HBV 模型 的 融雪 径流 模拟 效果 比 站 点 数据 
驱动 HBV 模型 的 融雪 径流 模拟 效果 更 好 , 格 点 数据 
在 研究 区 进行 融雪 洪水 模拟 的 适用 性 更 强 。 

结合 2 套数 据 的 对 比 结果 ,对 比 2 套 数据 驱动 
HBV 模型 模拟 出 的 融雪 洪水 过 程 与 实测 洪水 过 程 ， 
分 析 2 套 数据 在 研究 区 的 精度 与 误差 情况 。 在 每 年 
夏季 的 洪峰 部 分 ,2 套数 据 驱 动 HBV 模型 模拟 出 融 
雪 洪 水 的 洪峰 均 小 于 实测 洪水 的 洪峰 ,站 点 数据 的 
这 一 问题 更 为 突出 ,站 点 数据 模拟 出 的 这 里 的 洪峰 
流量 低 于 格 点 数据 模拟 出 的 流量 。 实 测 洪水 过 程 
中 ,洪峰 前 的 涨 洪 过 程 往往 比较 迅速 ,如 2008 年 7 
H 28—29 日 ,实测 流量 由 60.7 m-s KA 103 m-s”, 
2009 4E. 7 H 13—14 日 ,实测 流量 由 60.1 m? s" HRA 
111 m’+s',20114E6 H 18—19 日 ,实测 流量 由 57.4 
ms HEN 116 ms!。 而 2 套数 据 模 拟 出 的 融雪 洪 


水 流量 的 涨 洪 过 程 则 相对 平缓 。 对 比 相应 时 段 的 
气温 、 降 水 数据 ,2008 年 7 月 28 日 格 点 数据 降水 量 
为 11.4 mm, 站 点 数据 降水 量 则 仅 有 4.2 mm ; 2009 年 
7 月 13 日 格 点 数据 与 站 点 数据 的 降水 量 分 别 为 6.6 
mm 5 4.1 mm;2011 年 6 月 18 日 , 格 点 数据 降水 量 为 
21.0 mm, 站 点 数据 降水 量 为 10.0 mm。 根 据 模 型 结 
构 ,夏季 站 点 数据 驱动 HBV 模 型 模拟 出 的 洪峰 流量 
过 低 可 能 的 原因 是 ,在 夏季 ,所 用 站 点 的 降水 数据 
偏 低 或 者 所 用 站 点 的 气温 数据 偏 低 ,结合 2.1 的 分 
析 , 事 实 上 ,在 夏季 ,所 用 站 点 数据 的 气温 普遍 高 于 
格 点 数据 。 由 此 可 见 ,在 夏季 ,站 点 数据 驱动 HBV 
模型 模拟 出 的 洪峰 流量 过 低 的 原因 是 所 用 站 点 的 
降水 数据 偏 低 的 可 能 性 更 大 。 此 外 , 格 点 数据 驱动 
HBV 模型 模拟 出 的 洪峰 流量 也 低 于 实测 流量 ,只 有 
对 应 着 较 大 降水 的 情况 下 ,模拟 的 洪水 过 程 才 体 现 
出 陡 涨 的 过 程 ,因此 2 套 降水 数据 在 夏季 的 降水 量 
大 的 部 分 的 精度 均 有 待 提 高 ,站 点 数据 在 降水 量 大 
的 部 分 的 精度 比 格 点 数据 低 。 

2008 年 上 半年 ,2 套数 据 驱 动 HBV 模型 模拟 出 
融雪 洪水 的 流量 均 小 于 实测 值 ,主要 原因 是 模型 状 
态 变 量 的 初始 值 为 0, 需 要 运行 一 段 时 间 。 总 体 来 
看 ,除了 2009 年 ,2010 年 ,2011 年 的 春季 ,2009 年 的 
8—10 H ,2011 年 的 秋季 外 ,在 2008 一 2011 年 的 其 
他 时 段 ,用 站 点 数据 驱动 HBV 模 型 模拟 出 的 融雪 洪 
水 流量 普遍 小 于 实测 值 。 根 据 2.1 中 站 点 数据 与 格 
点 数据 对 比分 析 的 结果 ,推测 可 能 的 原因 是 所 用 站 
点 降水 数据 偏 低 , 导 致 模拟 出 的 融雪 洪水 流量 
偏 低 。 

在 2009 年 .2010 年 .2011 年 的 4 一 5 月 ,站 点 数 
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图 3 融合 数据 与 2 套数 据 驱 动 下 的 模拟 结果 


Fig.3 Simulation results driven by syncretic data and two sets of data 


据 驱 动 HBV 模型 模拟 出 的 融雪 洪水 流量 普遍 远大 
于 实测 流量 过 程 , 且 会 模拟 出 若干 洪峰 ,而 实测 流 
量 过 程 则 没有 这 些 洪峰 。 根 据 2 套 数据 对 比分 析 的 
结果 和 HBV 模 型 的 结构 ,推测 其 主要 原因 是 4 月 往 
往 处 于 一 年 中 气温 刚 开 始 持续 大 于 融雪 国 值 温度 
的 时 候 , 气 温 升 高 ,而 站 点 气温 数据 俩 高 ,因此 导致 
从 4 月 开始 模型 计算 出 的 融雪 量 人 高 ,从 而 计算 出 
的 春季 融雪 洪水 流量 偏 高 。 相 比 之 下 , 格 点 数据 驱 
动 HBV 模 型 模拟 出 的 洪水 流量 过 程 的 峰 现时 间 与 
实测 流量 过 程 总 体 较为 吻合 。 

在 每 年 的 冬季 ,实测 流量 往往 在 1.8~3.8 ms 
左右 , 格 点 数据 与 站 点 数据 驱动 HBV 模型 模拟 出 的 
冬季 流量 过 程 均 小 于 实测 流量 过 程 ,站 点 数据 驱动 
HBV 模 型 模拟 出 的 冬季 流量 更 小 ,基本 小 于 2 mss, 


格 点 数据 驱动 HBV 模型 模拟 出 的 冬季 流量 往往 维 
持 在 1.5~2.5 ms 左右 , 当 冬 季 实 测 流量 在 2 m? s" 
左右 时 , 格 点 数据 驱动 HBV 模 型 模拟 出 的 冬季 流量 
与 实测 流量 过 程 较 为 吻合 。 

从 数据 集 的 制作 过 程 分 析 误 差 原 因 , 格 点 数据 
是 基于 国家 气象 信息 中 心 基础 资料 专项 最 新 整编 
的 中 国 地 面 高 密度 台 站 的 降水 .气温 资料 ,利用 
ANUSPLIN 软件 的 薄 盘 样 条 法 (TPS, ThinPlate 
Spline) 进 行 空间 插值 ,生成 中 国 地 面 水 平分 辨 率 
0.5°x0.5° 的 日 值 降水 与 气温 格 点 数据 。 而 本 文 的 站 
点 数据 所 选取 的 站 点 为 研究 区 外 的 3 个 邻近 气象 站 
与 研究 区 内 1 个 水 文 站 。 在 山区 ,降水 受 地 形 影 响 
较 大 ,有 很 强 的 空间 变异 性 ,气温 也 具有 海拔 高 度 
引起 的 垂直 差异 。 这 些 站 点 的 气温 与 降水 数据 在 
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研究 区 代表 性 比较 差 ,精度 较 低 ,但 是 在 研究 区 再 
没有 别 的 站 点 。 因 此 ,在 后 续 的 研究 中 ,为 减 小 融 
雪 潜 水 模拟 的 误差 ,可 以 根据 海拔 高 程 . 坡 向 坡度 
等 因素 ,对 临近 站 点 数据 进行 修正 。 
2.3 融雪 洪水 特征 分 析 

由 图 3 可 以 看 出 ,在 研究 区 ,每 一 年 中 最 大 的 洪 
水 均 出 现在 夏季 , 且 往 往 涨 洪 迅 速 ,洪峰 陡峭 。 实 
测 日 平均 流量 常常 出 现 某 一 天 突然 陡 涨 的 现象 。 
选取 2000 4E. , 2001 4E. , 2008 年 为 例 ,2000 年 7 月 
15—16 日 7 月 17 一 18 H ,2001 4E. 7 H 29—30 日 和 
2008 年 7 月 28 一 29 日 实测 流量 增长 均 接 近 2 倍 。 
年 第 1 场 洪峰 所 在 的 日 平均 流量 超过 50 m^ s" 的 洪 
水 往往 出 现在 夏季 。 夏 季 是 全 年 洪水 流量 最 大 的 
季节 ,2000 一 2011 年 夏季 石门 站 实测 日 平均 流量 的 
多 年 平均 值 为 43.9 mx…s'。2000 一 2011 年 间 ,最 大 
的 日 平均 流量 为 2007 年 7 月 10 日 的 171 meso X} 
于 洪峰 部 分 的 日 平均 流量 ,2 套数 据 驱动 HBV 模型 
模拟 出 的 流量 往往 偏 小 。 

在 春季 ,日 平均 流量 由 冬季 的 1.8~3.8 mss" c 
右 开 始 缓慢 上 升 。 受 温度 升 高 . 积 雪 融化 与 降水 的 
影响 ,部 分 年 份 ,如 2000 年 .2002 年 .2003 年 .2006 
年 .2007 年 .2008 年 .2009 年 在 4 月 底 或 5 月 会 形成 
较 小 的 洪峰 ,洪峰 位 置 的 日 平均 流量 在 20~35 m^ s" 
左右 。2000 一 2011 年 春季 实测 日 平均 流量 的 多 年 
平均 值 为 5.66 m «s! 。 站 点 数据 驱动 HBV 模型 模拟 
出 的 春季 流量 过 程 常 常 偏 大 。 

每 年 的 9 月 开始 ,秋季 位 于 一 年 中 最 后 1 场 较 
大 的 洪水 的 退 水 过 程 , 退 水 过 程 往往 比较 缓慢 ,部 
分 年 份 的 日 平均 流量 在 这 一 段 退 水 过 程 中 有 上 下 
波动 过 程 ,如 2007 年 9 月 上 旬 ,2009 年 9 月 上 旬 、 
20104E 9 H ER], 2000—2011 年 秋季 石门 站 实测 
日 平均 流量 的 多 年 平均 值 为 10.5 m-s 

冬季 流量 普遍 较 小 ,一 般 保持 在 1.8~3.8 ms 
左右 ,波动 很 小 。2000 一 2011 年 冬季 实测 日 平均 流 
量 的 多 年 平均 值 为 2.80 m-s" 


3 结论 

(1) 格 点 数据 驱动 下 水 文 模型 的 融雪 径流 模拟 
效果 总 体 优 于 利用 站 点 观测 数据 驱动 的 水 文 模型 
的 模拟 效果 。 格 点 数据 驱动 下 水 文 模型 模拟 融雪 
径流 的 Nash 系数 在 验证 期 为 0.792, 高 于 站 点 数据 
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驱动 下 的 0.433。 格 点 数据 驱动 下 水 文 模型 模拟 融 
雪 径流 的 平均 相对 误差 为 2.44%, 低 于 站 点 数据 驱 
动 下 的 14.64%。 基 于 国家 级 高 密度 台 站 插值 而 来 
的 格 点 数据 在 研究 区 融雪 洪水 模拟 的 适用 性 比 本 
文 所 选取 的 研究 区 外 的 3 个 邻近 气象 站 与 研究 区 内 
1 个 水 文 站 的 气温 与 降水 数据 更 高 。 

(2) 分 析 模 拟 与 实测 的 融雪 洪水 特征 ,在 大 潜 
水 的 部 分 ,实测 洪水 过 程 涨 洪 迅 速 ,2 套 数据 模拟 出 
的 洪峰 流量 都 普遍 偏 小 ,站 点 数据 驱动 HBV 模 型 模 
拟 出 的 融雪 洪水 洪峰 部 分 的 误差 大 于 格 点 数据 驱 
动 HBV 模 型 的 模拟 结果 。 在 春季 ,站 点 数据 驱动 
HBV 模 型 模拟 出 的 流量 过 程 比 实测 过 程 偏 大 ,这 主 
要 是 因为 所 用 站 点 的 气温 数据 偏 高 。 在 每 年 的 11 
月 下 旬 至 次 年 2 月 上 名 ,2 套数 据 驱 动 HBV 模型 得 
到 的 流量 模拟 值 都 普遍 低 于 实测 值 。 对 于 格 点 数 
据 , 当 实测 流量 在 2 ms 左右 时 ,模拟 值 与 实测 值 
会 比较 吻合 。 冬 季 站 点 数据 在 研究 区 的 精度 则 有 
待 提高 。 

(3) 结合 2 套数 据 的 对 比 结果 与 融雪 径流 的 模 
拟 情况 可 知 ,本 文 所 选取 的 研究 区 外 的 3 个 邻近 气 
象 站 与 研究 区 内 1 个 水 文 站 的 气温 与 降水 数据 在 研 
究 区 的 精度 较 低 ,特别 是 日 降水 数据 的 算术 平均 值 
在 夏季 偏 低 , 日 平均 气温 算术 平均 值 在 春季 偏 高 。 
因为 在 山区 ,降水 受 地 形 影 响 较 大 ,有 很 强 的 空间 
变异 性 ,气温 也 具有 海拔 高 度 引 起 的 垂直 差异 。 这 
些 站 点 的 气温 与 降水 数据 在 研究 区 代表 性 比较 差 ， 
精度 较 低 ,但 在 人 研究 区 无 其 他 站 点 。 本 文 仅 基于 水 
文 模型 模拟 出 的 融雪 洪水 径流 与 实测 径流 对 比 , 分 
别 对 2 套数 据 集 进行 了 对 比分 析 及 其 融雪 型 径流 模 
拟 的 适用 性 分 析 。 在 后 续 的 研究 中 ,可 以 根据 海拔 
高 程 . 坡 向 坡度 等 因素 ,对 临近 站 点 数据 进行 修正 ， 
以 减 小 融雪 洪水 模拟 的 误差 。 
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Applicability assessment of two meteorological datasets in areas lacking data 
with the Hutubi River Basin as an example 
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Abstract: In ungauged areas, a lack of observational data (e.g., precipitation and temperature data) can affect the 
accuracy of flood forecasting, which in turn affects the implementation of flood control, drought relief, and water 
resources planning. Therefore, it is necessary to conduct applicability analysis of existing precipitation and 
temperature data in ungauged areas. By applying the HBV (Hydrologiska Byrans Vattenbalansavdelning) model, 
including a snowmelt runoff simulation, to the area upstream of the Shimen hydrological station in the Hutubi 
River Basin, and based on a comparison of simulated snowmelt flood runoff and measured runoff, a dataset 
including gridded daily precipitation and temperature in China and observational data from meteorological 
stations were analyzed and compared. The applicability of the snowmelt runoff simulation was also analyzed. 
Results showed that the snowmelt runoff simulation was more effective when used with the hydrological model 
driven by grid data than with the hydrological model driven by station observation data. The Nash coefficients of 
the snowmelt runoff from the hydrological model driven by grid data and the hydrological model driven by 
station observation data were 0.792 and 0.433, respectively, within the verification period. In addition, some 
characteristics of the snowmelt floods in the Hutubi River Basin were analyzed and the possible causes of errors 
in the simulated snowmelt runoff were also assessed based on differences in driving data. These results can help 
improve the accuracy of snowmelt flood forecasts in areas that lack data on snowmelt characteristics. 

Keywords: precipitation data; ungauged area; snowmelt flood; runoff forecast; applicability assessment; Hutu- 
bi River Basin 


